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Abstract 



Feedback loops play an important role in determining the dynamics of bio- 
logical networks. In order to study the role of negative feedback loops, this 
paper introduces the notion of "distance to positive feedback (PF-distance)" 
which in essence captures the number of "independent" negative feedback 
loops in the network, a property inherent in the network topology. Through 
a computational study using Boolean networks it is shown that PF-distance 
has a strong influence on network dynamics and correlates very well with 
the number and length of limit cycles in the phase space of the network. To 
be precise, it is shown that, as the number of independent negative feed- 
back loops increases, the number (length) of limit cycles tends to decrease 
(increase). These conclusions are consistent with the fact that certain nat- 
ural biological networks exhibit generally regular behavior and have fewer 
negative feedback loops than randomized networks with the same numbers 
of nodes and connectivity. 

Key words: Feedback loops; Boolean networks; gene regulatory net- 
works; limit cycles 
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Introduction 

An understanding of the design principles of biochemical networks, such as 
gene regulatory, metabolic, or intracellular signaling networks is a central 
concern of systems biology. In particular, the intricate interplay between 
network topology and resulting dynamics is crucial to our understanding of 
such networks, as is their presumed modular structure. Features that relate 
network topology to dynamics may be considered "robust" in the sense that 
their influence does not depend on detailed quantitative features such as 
exact flux rates. A topological feature of central interest in this context 
is the existence of positive and negative feedback loops. There is broad 
consensus that feedback loops have a decisive effect on dynamics, which 
has been studied extensively through the analysis of mathematical network 
models, both continuous and discrete. Indeed, it has long been appreciated 
by biologists that positive and negative feedback loops play a central role in 
controlling the dynamics of a wide range of biological systems. Thomas et 
al. Thomas et al. (1) conjectured that positive feedback loops are necessary 
for multistationarity whereas negative feedback loops are necessary for the 
existence of periodic behaviors. Proofs for different partial cases of these 
conjectures have been given, see, e.g., Soule (j^), Plahte et al. (3), Gauze 
(jj), Cinquin and Demongeot (jsl), Angeli et al. Moreover, it is widely 
believed Sontag (jTI) that an abundance of negative loops should result in the 
existence of "chaotic" behavior in the network. This paper provides strong 
evidence in support of this latter conjecture. 

We focus here on Boolean network (BN) models, a popular model type 
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for biochemical networks, initially introduced by S. Kauffman KaufFman (|8|). 
In particular, we study BN models in which each directed edge can be char- 
acterized as either an inhibition or an activation. In Boolean models of bio- 
logical networks, each variable can only attain two values (0/1 or "on/off"). 
These values represent whether a gene is being expressed, or the concen- 
tration of a protein is above a certain threshold, at time t. When detailed 
information on kinetic rates of protein-DNA or protein-protein interactions 
is lacking, and especially if regulatory relationships are strongly sigmoidal, 
such models are useful in theoretical analysis, because they serve to focus 
attention on the basic dynamical characteristics while ignoring specifics of 



reaction mechanisms, see Kauffman (|9|), Kauffman and Glass did), Albert 
and Othmer (11), Chaves et al. (Q)). 

Boolean networks constructed from monotone Boolean functions (i.e. 
each node or "gate" computes a function which is increasing on all argu- 
ments) are of particular interest, and have been studied extensively, in the 



electronic circuit design and pattern recognition literature Gilbert (|l3l ) , Min 



sky 



as well as in the computer science literature; see e.g. Aracena et al. 
id). Goles and Hernandez jl?! ) for recent references. For Boolean and 



all other finite iterated systems, all trajectories must either settle into equi- 
libria or to periodic orbits, whether the system is made up of monotone 
functions or not, but monotone networks have always somewhat shorter cy- 
cles. This is because periodic orbits must be anti-chains, i.e. no two different 
states can be compared; see Gilbert jisl ). Smith (18). An upper bound may 



be obtained by appealing to Sperner's Theorem (Anderson (jl9l )): Boolean 



systems on n variables can have orbits of period up to 2", but monotone 
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systems cannot have orbits of size larger than ( ^JJ^^ ) 2" y ^ ; these are 
all classical facts in Boolean circuit design Gilbert (|l3l). It is also known 
that the upper bound is tight Gilbert (jl3i ). in the sense that it is possible to 
construct Boolean systems on n variables, made up of monotone functions, 
for which orbits of the maximal size (|^„"2j) Si'^^n by Sperner's Theorem 
exist. This number is still exponential in n. However, anecdotal experi- 
ence suggests that monotone systems constructed according to reasonable 
interconnection topologies and/or using restricted classes of gate functions, 
tend to exhibit shorter orbits Tosic and Agha Greil and Drossel (12 ll ). 



One may ask if the architecture of the network, that is, the structure of 
its dependency (also called interconnection) graph, helps insure shorter or- 



bits. In this direction, the paper Aracena et al. (|l5l ) showed that on certain 



graphs, called there "caterpillars" , monotone networks can only have cycles 
of length at most two in their phase space. 

The present paper asks the even more general question of whether net- 
works that are not necessarily made up from monotone functions, but which 
are "close to monotone" (in a sense to be made precise, roughly meaning that 
there are few independent negative loops) have shorter cycles than networks 

which are relatively farther to monotone. 

] 

In Sontag Cil), we conjectured that "smaller distance to monotone" should 
correlate with more ordered (less "chaotic") behavior, for random Boolean 
networks. A partial confirmation of this conjecture was provided in Kwon 
and Cho (22), where the relationship between the dynamics of random 
Boolean networks and the ratio of negative to positive feedback loops was 
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investigated, albeit only for the special case of small Kauffman-type NK 
and NE networks, and with the additional restriction that all nodes have 
the same function chosen from AND, OR, or UNBIAS. Based on computer 



simulations, the authors of Kwon and Cho (|22l ) found a positive (negative) 
correlation between the ratio of fixed points (other limit cycles) and the ratio 
of positive feedback loops. Observe that this differs from our conjecture in 
two fundamental ways: (1) our measure of disorder is related to the number 
of "independent" negative loops, rather than their absolute number, and 
(2) we do not consider that the number of positive loops should be part of 
this measure: a large number of negative loops will tend to produce large 
periodic orbits, even if the negative to positive ratio is small due to a larger 
number of positive loops. 

Thus, in the spirit of the conjecture in Sontag (j3), the current paper 
has as its goal an experimental study of the effect of independent negative 
feedback loops on network dynamics, based on an appropriately defined 
measure of distance to positive-feedback. We study the effect of this distance 
on features of the network dynamics, namely the number and length of limit 
cycles. Rather than focusing on the number of negative feedback loops in the 
network, as the characteristic feature of a network, we focus on the number 
of switches of the activation/inhibition character of edges that need to be 
made in order to obtain a network that has only positive feedback loops. We 
relate this measure to the cycle structure of the phase space of the network. 
It is worth emphasizing that the absolute number of negative feedback loops 
and the distance to positive feedback are not correlated in any direct way, 
as it is easy to construct networks with a fixed distance to positive feedback 
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that have arbitrarily many negative feedback loops, see Figure [H 
Motivations 

There are three different motivations for posing the question that we ask in 
this paper. The first is that most biological networks appear to have highly 
regular dynamical behavior, settling upon simple periodic orbits or steady 
states. The second motivation is that it appears that real biological networks 



such as gene regulatory networks anc 
close to monotone DasGupta et al. 



OTotein signaling networks are indeed 



23|), Sontag ([2^), Maayan et al. 



Thus, one may ask if being close to monotone correlates in some way with 
shorter cycles. Unfortunately, as mentioned above, one can build networks 
that are monotone yet exhibit exponentially long orbits. This suggests that 
one way to formulate the problem is through a statistical exploration of 
graph topologies, and that is what we do here. 

A third motivation arises from the study of systems with continuous 
variables, which arguably provide more accurate models of biochemical net- 
works. There is a rich theory of continuous-variable monotone (to be more 
precise, "cooperative") systems. These are systems defined by the prop- 
erty that an inequality a(0) < b(0) in initial conditions propagates in 
time so that the inequality a(t) < b(t) remains true for all future times 
t > 0. Note that this is entirely analogous to the Boolean case, when one 
makes the obvious definition that two Boolean vectors satisfy the inequality 
a = (oi, . . . , an) < b = . . . , 6„) if < bi for each i = 1, . . . , n (setting 
< 1). Monotone continuous systems have convergent behavior. For ex- 
ample, in continuous-time (ordinary differential models), they cannot admit 
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any possible stable oscillations Hadeler and Glas (26), Hirsch 



Hirsch 



and Smith (|28l). and, when there is only one steady state, every bounded 



solution converges to this unique steady state (monostability) , see Dancer 



Dancer (|29l ). When, instead, there are multiple steady-states, the Hirsch 



Generic Convergence Theorem Smith (llSl ). Hirsch and Smith (|28l ). Hirsch 



l3ll ) is the fundamental result; it states, under an additional technical 
assumption ("strong" monotonicity) that generic bounded solutions must 
converge to the set of steady states. For discrete-time strongly monotone sys- 
tems, generically also stable oscillations are allowed besides convergence to 
equilibria, but no more complicated behavior. In neither case, discrete-time 
or continuous-time continuous monotone systems, one observes "chaotic" 
behavior. It is an open question whether continuous systems that are in 
some sense close to being monotone have more regular behavior, in a sta- 
tistical sense, than systems that are far from being monotone, just as for 
the Boolean analog considered in this paper. The Boolean case is more 
amenable to computational exploration than continuous- variable systems, 
however. Since long orbits in discrete systems may be viewed as an analog 
of chaotic behavior, we focus on lengths of orbits. 

One can proceed in several ways to define precisely the meaning of dis- 
tance to positive feedback. One associates to a network made of unate (defi- 
nition below) gate functions a signed graph whose edges have signs (positive 
or negative) that indicate how each variable affects each other variable (acti- 
vation or inhibition). The first definition, explored in Sonta g (17 ), DasGupta 



et al. (1231 ), Sontag (j24l ). Dasgupta et al. (|32| ). Hiiffner et al. (1331 ) starts from 



the observation that in a network with all monotone node functions there are 
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no negative undirected cycles. Conversely, if the dependency graph has no 
undirected negative parity cycles (a "sign-consistent" graph), then a change 
of coordinates (globally replacing a subset of the variables by their com- 
plements) renders the overall system monotone. Thus, asking what is the 
smallest number of sign-flips needed to render a graph sign-consistent is one 
way to define distance to monotone. This approach makes contact with 
areas of statistical physics (the number in question amounts to the ground 
energy of an associated Ising spin-glass model), as well as with the general 



theory of graph-balancing for signed graphs Zaslavsky (1341 ) that originated 



with Harary Harary (|35l). It is also consistent with the generally accepted 
meaning of "monotone with respect to some orthant order" in the ODE 
literature as a system that is cooperative under some inversion of variables. 

A second, and different, definition, starts from the fact that a network 
with all monotone node functions has, in particular, no negative-sign directed 
loops. For a strongly connected graph, the property that no directed neg- 
ative cycles exist is equivalent to the property that no undirected negative 
cycles exist. However, for non-strongly connected graphs, the properties are 
not the same. Thus, this second property is weaker. The second property is 
closer to what biologists and engineers mean by "not having negative feed- 
backs" in a system, and hence is perhaps more natural for applications. In 
addition, it is intuitively clear that negative feedbacks should be correlated 
to possible oscillatory behavior. (This is basically Thomas' conjecture. See 



Sontag (l2J) for precise statements for continuous-time systems; interestingly, 
published proofs of Thomas' conjecture use the first definition, because they 
appeal to results from monotone dynamical systems.) Thus, one could also 
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define distance to monotone as the smallest number of sign-flips needed to 
render a graph free of negative directed loops. To avoid confusion, we will 
call this notion, which is the one studied in this paper, distance to positive- 
feedback, or just "PF-distance" . 

Theory 

Distance to positive-feedback 

We give here the basic definitions of the concepts relevant to the study. 

Definition 1 Let k = {0, 1} be the field with two elements. We order the 
two elements as < 1. This ordering can be extended to a partial ordering 
on A;" by comparing vectors coordinate-wise in the lexicographic ordering. 

1. A Boolean function h : /c" — > k is monotone if, whenever a < b 
coordinate-wise, for a, b G k^, then h{a.) < h{h). 

2. A Boolean function h is unate if, whenever Xi appears in h, the fol- 
lowing holds: Either 

(a) For all ai, . . . ,aj_i,ai+i, . . . ,a„ G A;, 



h{ai, . . . ,ai_i,0,aj+i, ...,an) < 



h{ai, . . . ,ai-i, l,ai+i, . 



,an), or 



(b) For allai,...,ai-i,ai+i,..., 



Qn & k 



h{ai, . . . , Oj-i, 0, aj+i, . . . 



h{ai, . . . , tti-i, 1, Oj+i, . 
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The definition of unate function is equivalent to requiring that whenever 
ai appears in h, then it appears either everywhere as ai or everywhere as 
-^Qi := 1 + aj. 

Let / be a Boolean network with variables xi, . . . and coordinate 
functions /i, . . . , /„. That is, / = (/i, . . . , /„) : A;" — > /c". We can associate 
to / its dependency graph T>{f): The vertices are fi, . . . , corresponding to 
the variables xi, . . . , x„, and there is an edge Vi — > Vj if and only if Xi appears 
in fj. If all coordinate functions fi of / are unate, then the dependency graph 
of / is a signed graph. Namely, we associate to an edge Vi Vj a "+" if fj 
preserves the ordering as in 2(a) of Definition [1] and a "-" if it reverses the 
ordering as in 2(b) of Definition [TJ For later use we observe that this graph 
(as any directed graph) can be decomposed into a collection of strongly 
connected components, with edges between strongly connected components 
going one way but not the other. (Recall that a strongly connected directed 
graph is one in which any two vertices are connected by a directed path.) 
That is, the graph can be represented by a partially ordered set in which the 
strongly connected components make up the elements and the edge direction 
between components determines the order in the partially ordered set. 

Definition 2 Let f be a Boolean network with unate Boolean functions and 
T){f) he its signed dependency graph. Then 

1. f is a positive-feedback network (PF) if T>{f) does not contain any 
odd parity directed cycles. ( The parity of a directed cycle is the product 
of the signs of all the edges in the cycle.) 
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2. The PF-distance of f is the smallest number of signs that need to be 
changed in the dependency graph to obtain a PF network. We denote 
this number by \'D{f)\ or simply \ f\. 

Notice that for a given directed graph G, different assignments of sign to 
the edges produce graphs with varying PF-distance. In particular, there is 
a maximal PF-distance that a given graph topology can support. 

The dynamics of / are presented in a directed graph, called the phase 
space of /, which has the 2" elements of fc" as a vertex set, and there is an 
edge a ^ b if /(a) = b. It is straightforward to see that each component 
of the phase space has the structure of a directed cycle, a limit cycle., with 
a directed tree feeding into each node of the limit cycle. The elements of 
these trees are called transient states. 

In this paper we relate the dynamics of a Boolean network to its PF- 
distance. The following is a motivational example that explains the main 
results. 

Example 3 Let G he the directed graph depicted in Fig. (left). It is easy 
to check that the maximal PF-distance of G is 3. Let f = (xa V -'X4, xi A 
X2, X2 A -ia;4, ^x^) : {0, 1}^ — > {0, 1}^ and g = {^x^ V x^, xi A ^X2,X2 A 
3^4) "■2^3) '■ {0, 1}'^ — > {0, 1}^. It is clear that f and g are sign-modifications 
of the same PF network (xa V X4,xi A X2,X2 A X4,X3), in particular, they 
have the same (unsigned) dependency graph. However, the PF-distance of f 
is while it is 3 for g. The phase space of f is depicted in Fig. (middle) 
and that of g is on the right. Notice that f has two limit cycles of lengths 1 
and 2, respectively, while g has only one limit cycle of length 4. 
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For each distance < d < 3, we analyze the dynamics of 10 random 
PF networks and their sign modifications of distance d on the directed graph 
in Fig. (left). The average of the numbers (lengths) of limit cycles is 
computed as in Table [IJ The best fit-line of the averages of the number 
(length) of limit cycles is computed and its slope is reported as in Fig. O 
The details of this analysis are provided in the Supporting Inforamtion. 

We have repeated the experiment in Example [3] above for many different 
graphs and observed that the slope of the best fit-hne of the length (resp. 
number) of hmit cycles is positive (resp. negative) most of the time. In the 
methods section we present the details of the experiment and the algorithms 
used in the computations. The results of these experiments are described 
next. 

Methods 

The main results of this paper relate the PF-distance of Boolean networks 
with the number and length of their limit cycles. Specifically, our hypoth- 
esis is that, for Boolean networks consisting of unate functions, as the PF- 
distance increases, the total number of limit cycles decreases on average and 
their average length increases. This is equivalent to saying that for most or 
all experiments this slope is negative for the number of limit cycles and is 
positive for their length. 

To test this hypothesis we analyzed the dynamics of more than six million 
Boolean networks arranged in about 130,000 experiments on random graphs 
with the number of nodes 5, 7, 10, 15, 20, or 100 and maximum in-degree 5 



Dynamics of Boolean networks 



13 



for each node. 

Random generation of unate functions 

We generated a total of more than 130,000 random directed graphs, where 
each graph has 5,7,10,15, 20, or 100 nodes, with maximum in-degree 5 for 
each node. The graphs were generated as random adjacency matrices, with 
the restriction that each row has at least one 1 and at most five I's. For each 
graph directed G, we generated 10 Boolean networks with unate functions 
and dependency graph G, by using the following fact. 

Lemma 4 A Boolean function f of n variables is unate if and only if it is 
of the form /(x) = gi(x + s), where g is a monotone Boolean function of n 
variables and s G A;" and "+ " denotes addition modulo 2. 

Proof. If / is unate then each variable Xj appears in / always as xi or 
always as -iXj. Suppose that all Xi appear without negations. Then / is 
constructed using A and V. Hence / is monotone. Otherwise, let s € A;" be 
the vector whose ith entry is 1 if and only if xi appears as -iXj in /. Then 
^(x) = /(x + s) is a monotone function and /(x) = ^(x + s). The converse 
is clear. 

So in order to generate unate functions it is sufficient to generate mono- 
tone functions. We generated the set Mj of monotone functions in i variables 
by exhaustive search for z = 1, . . . , 5. (For example, M5 has 6894 elements.) 
Unate functions for a given signed dependency graph can then be generated 
by choosing random functions from Mj and random vectors s G fc". The 
nonzero entries in s for a given node correspond to the incoming edges with 
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negative sign in the dependency graph. Using this process we generated 
Boolean networks with unate Boolean functions. 
We then carry out the following steps. 

The Experiment 

Let G be a random unsigned directed graph on n nodes with a maximal 
PF-distancc t, and let D < t. Consider 10 unate Boolean networks chosen 
at random with G as their dependency graph. 

1. For 1 < d < D, let Gd be a signed graph of G of distance d. 

(a) For each network / of the ten networks, 

i. Let g he a modified network of / such that 'D{g) = Gd', the 
signed dependency graph of ^ is G^. 

ii. Compute the number and length of all limit cycles in the 
phase space of g. 

(b) Compute the average number N (resp. average length L) of limit 
cycles in the phase spaces of the g's. 

2. Compute the slope sn (resp. sl) of the best fit-line of the N's (resp. 
L's). 



Computation of PF-distance 

Let / be a Boolean network with unate Boolean functions and let |/| be its 
PF-distance. The proofs of the following facts are straightforward. 
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1. Suppose the dependency graph of / has a negative feedback loop at a 
vertex. Let /' be the Boolean network obtained by changing a single 
sign to make the loop positive. Then |/| = |/'| + 1. 

2. Let Hi, . . . Hg be the strongly connected components of the depen- 
dency graph D{f). Then 



\V{f)\ = Y,\H, 



i=l 



The algorithm for computing now follows. 



Algorithm: Distance to PF 

Input: A signed, directed graph G. 
Output: |G|; the PF-distance of G. 

Let d = 0. 

1. Let Gi,...,Gr be the collection of all 
signed graphs obtained by making ex- 
actly d sign changes in G. 



2. For i = 1, . . . ,r 



If Gi is PF, then RETURN |G| = d. 

3. Otherwise, d := d + 1, Go to Step ([T]) 
above. 



In Step ([2]) above, to check whether a strongly connected graph is PF, 
it is equivalent to check whether it has any (undirected) negative cycles. 



which can easily be done in many different ways, see, e.g., Sontag (j24l ). This 
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algorithm must terminate, since G has finitely many edges and hence the 
PF-distance of G is finite. 

If G has m directed edges, then there are 2™ possible sign assignments. 
However, to compute the maximal PF-distance, one does not need to find 
the PF-distance of such possible assignments, see Supporting Information 
for the algorithm we used to compute the maximal PF-distance. 

Results 

In Table [2] we present the percentage of experiments that do conform to our 
hypothesis for the average of number of limit cycles as well as for the average 
length of limit cycles. It should be mentioned here that a computationally 
expensive part of an experiment is the computation of the maximal PF- 
distance of a given directed graph and becomes prohibitive for even modest- 
size graphs, with, e.g., 10 nodes. So unlike in Example [3l for networks on 
more than 5 nodes, we only considered PF-distances that are less than or 
equal to the number of nodes in the network. (See the Methods Section for 
a detailed description of the experiment.) In fact for graphs with 20 (resp. 
100) nodes, all considered networks have PF-distance less than or equal to 
5 (resp. 10). We argue below that this is the reason for the drop in the 
percentage of experiment that conform to our hypothesis as the number of 
nodes increases. 

For networks on 5 nodes, we analyzed the dynamics of 4000 experiments 
by varying the PF-distance considered in the computations. Table [3] shows 
the number of experiments that do not conform to our hypothesis as we 
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vary the considered PF-distance. We also present the results in the form 
of histograms, where the horizontal axis represents the slope of the hnes of 
best fit and the vertical axis represents the percentage of experiments that 
confirm our hypothesis. Figure H] shows the results of the 4000 experiments 
on 5-node networks. The histograms, from left to right, show the results 
when PF-distance of the network is 25%, resp. 50%, resp. 75%, resp. 
100% of the maximal PF-distance. It can be seen in the right-most figure 
that almost all experiments show positive slope of the best-fit line, thereby 
conforming to the conjecture. Similar results for the average number of 
limit cycles are shown in Figure O demonstrating that if the PF-distance 
of networks is allowed the whole possible range, almost all the experiments 
conform to our hypothesis as we already noticed in Table O However, these 
histograms show the distribution of slopes over different distances. 

We have carried out similar computations for networks with 7 (5000 ex- 
periments)and 10 (6000 experiments) nodes; see the Supporting information 
for details. The results there are not quite as clear as for 5 node networks. 
For instance, for networks with 10 nodes (and up to 4 incoming edges per 
node) 31 out of 1000 experiments did not conform to our hypothesis for 
networks with PF-distance up to 5. 

In summary, the exhaustive computations confirm our hypothesis that, 
as the PF-distance increases, the total number of limit cycles decreases on av- 
erage and their average length increases. Furthermore, the slopes of the best- 
fit lines increasingly conform to our hypothesis the closer the PF-distance of 
the networks comes to the maximum PF-distance of the network topology. 
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Supporting Information 

For details of our analysis of the 5-node networks and all other considered 
networks, see Supporting Text. 



Discussion 

Negative feedback loops in biological networks play a crucial role in control- 
ling network dynamics. The new measure of "distance to positive feedback 
(PF-distance)" introduced in this paper is designed to capture the notion of 
"independent" feedback loops. We have shown that PF-distance correlates 
very well with the average number and length of limit cycles in networks, key 
measures of network dynamics. By analyzing the dynamics of more than six 
millions Boolean networks, we have provided evidence that networks with a 
larger number of independent negative feedback loops tend to have longer 
limit cycles and thus may exhibit more "random" or "chaotic" behavior. 
Furthermore, the number of limit cycles tends to decrease as the number of 
independent negative feedback loops increases. 

In general, the problem of computing the PF-distance of a network is 
NP-complete, as MAX-CUT can be mapped into it as a special case; see 



DasGupta et al. (1231 ). Sontag (24), Dasgupta et al. (|32l ) for a discussion for 
the analogous problem of distance to monotone. The question of computing 
distance to monotone has been the subject of a few recent papers DasGupta 



et al. (1231 ) . Dasgupta et al. (|32l ). Hiiffner et al. (|33l ). The first two of these 
proposed a randomized algorithm based on a semi-definite programming 
relaxation, while the last one suggested an efficient deterministic algorithm 
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for graphs with small distance to monotone. Since a strongly connected 
component of a graph is monotone if and only if it has the PF property, 
methods for computing PF distance for large graphs may be developed by 
similar techniques. Work along these lines is in progress. 
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Figure Legends 

Figure [1]. 

A graph that has an arbitrary number of negative loops, as many as the 
number of nodes in the second layer, but PF distance one: to avoid negative 
feedback, it suffices to switch the sign of the single (negative) arrow from 
the bottom to the top node. All unlabeled arrows are positive. 

Figure [2]. 

The dependency graph (left), the phase space of / (middle) and the phase 
space of g (right) from Example [3l These graph were generated using DVD 



Jarrah et al. (|36l). 



Table [H 

The average of the numbers (lengths) of limit cycles of the networks from 
Example [3l 

Table El 

The percentage of experiments that conform the hypotheses 
Table 

The number of experiments that did not conform to our hypothesis for 5- 
node networks. We considered PF-distance 25%, 50%, 75%, and 100% of the 
maximal distance . For each d, we considered 1000 experiments. 
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Figure [U. 

5-node networks. Histogram of slopes of best-fit lines to the average length 
of limit cycles (horizontal axis) vs. percentage of experiments with a given 
slope (vertical axis). The panels from left to right include networks with 
increasing PF-distance, with 25%, 50%, 75%, and 100% of the maximal 
distance. 

Figure [5]. 

5-node networks. Histogram of slopes of best-fit lines to the average number 
of limit cycles (horizontal axis) vs. percentage of experiments with a given 
slope (vertical axis). The panels from left to right include networks with 
increasing PF-distance, with 25%, 50%, 75%, and 100% of the maximal 
distance. 

Figure [6]. 

Histogram of slopes of best-fit lines to the average number, resp. length, 
of limit cycles (horizontal axis) vs. percentage of experiments with a given 
slope (vertical axis). The left two panels show the results for 15-node net- 
works, the other two panels those for 20-node networks. 
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Table 1: 



d 


Av. Num. 


Av. Len. 





3.5 


1.23 


1 


2.80 


1.25 


2 


2.50 


1.52 


3 


1.20 


3.50 
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amies of Boolean networks 



Table 2: 



n 


Num. of Exp. 


Av. Num. 


Av. Len. 


5 


117000 


99.75 


99.83 


7 


5000 


97.82 


99.92 


10 


6000 


95.70 


99.58 


15 


2921 


95.72 


98.25 


20 


331 


90.03 


94.86 


100 


659 


77.39 


93.93 



amies of Boolean networks 



Table 3: 



D 


Av. Num. 


Av. Len. 


Med. Num. 


Med. Len. 


25% 


26 


114 


29 


542 


50% 


4 


16 


18 


59 


75% 





1 


6 


3 


100% 


1 





2 
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Figure 5: 
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